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Abstract. We present a conjecture for the power-law exponent in the asymp- 
totic number of types of plane curves as the number of self-intersections goes 
to infinity. In view of the description of prime alternating links as fiype equiv- 
alence classes of plane curves, a similar conjecture is made for the asymptotic 
number of prime alternating knots. 

The rationale leading to these conjectures is given by quantum field theory. 
Plane curves are viewed as configurations of loops on a random planar lattices, 
that are in turn interpreted as a model of 2d quantum gravity with matter. 
The identification of the universality class of this model yields the conjecture. 

Since approximate counting or sampling planar curves with more than a few 
dozens of intersections is an open problem, direct confrontation with numeri- 
cal data yields no convincing indication on the correctness of our conjectures. 
However, our physical approach yields a more general conjecture about con- 
nected systems of curves. We take advantage of this to design an original and 
feasible numerical test, based on recent perfect samplers for large planar maps. 
The numerical data strongly support our identification with a conformal field 
theory recently described by Read and Saleur. 



1. Introduction. 

Our motivation for this work is the enumeration of topological equivalence classes 
of smooth open and closed curves in the plane (see Figure ^ precise definitions are 
given in Section 12.1(1 . The problem of characterizing closed curves was considered 
already by Gauss and has generated many works since then: see and references 
therein. Our interest here is in the numbers Op and ap of such open and closed 
curves with p self-intersections, and more precisely we shall consider the asymptotic 
properties of Op and ap when p goes to infinity. The numbers Op were given up to 
p = 10 in [HI and have been recently computed up to p = 22 by transfer matrix 
methods |1[)| . Asymptotically, as p goes to infinity, one expects the relation Up ^ 
Aap to hold (see below), so that we concentrate on the numbers Qp. 

In the present paper we propose a physical reinterpretation of the numbers Up 
that leads to the following conjecture, and we present numerical results supporting 
it. 

Conjecture 1. There exist constants t and c such that 

ap -ap CT^ -p^ , 

p-^oo 4 P — ''OO 

where 

(1) ^^_ l + \/l3 ^ _0 7g759_^ 
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Figure 1. An open plane curve and the associated closed curve. 



From the data of ^01 one has the numerical estimate: r = 11.4. But the main 
point in Conjectureniies not so much in the existence of r as in the explicit value of 
7. It should indeed be observed that 7, rather than t, is the interesting information 
in the asymptotic of a^. Observe for instance that the value of 7 is left unchanged 
if one redefines the size of a closed curve as the number p' = 2p of arcs between 
crossings. More generally, as discussed in Section|Sl the exponent 7 determines the 
branching behavior of generic large curves under the uniform distribution, and is 
universal in the sense that the same value is expected in the asymptotic of related 
families of objects like prime self-intersecting curves or alternating knots. 

Conjecturenis similar in nature to the conjecture of Di Francesco, Golinelli and 
Quitter on the asymptotic behavior of the number of plane meanders [Hj. The two 
problems do not fall into the same universality class (in particular the predictions 
for the exponent 7 are different in the two problems). However our approach to 
design a numerical test is applicable also to the meander problem. 

The rest of the paper is organized as follows. Precise definitions are given, 
and a more general family of drawings is introduced that play an important role 
in the identification of the associated physical model (Section EJ. The physical 
background leading to the conjecture is then discussed (Section|3I) and a numerically 
testable quantity is proposed (Section^]). The sampling method is briefly presented 
(Section 15)) before the analysis of numerical data (Section 01 . We conclude with 
some variants and corollaries of the conjecture (Section IHJ. 

2. Plane curves and colored planar maps 

2.1. Plane curves and doodles. For p a positive integer, let Ap be the set of 

equivalence classes of self-intersecting loops 7 in the plane, that is: (i) 7 is a smooth 
mapping R^; (ii) there are p points of self-intersection, all of which are 

regular crossings; (iii) two loops 7 and 7' are equivalent if there exists an orientation 
preserving homeomorphism h of the plane such that ^'{S^) = h{j{S^)). 

Similarly let Ap be the set of equivalence classes of self-intersecting open curves 
7 in the plane: (i) 7 is a smooth mapping [0, 1] and 7(0) and 7(1) belong to 

the infinite component of \ 7((0, 1)); (ii) there are p points of self-intersection, 
all of which are regular crossings; (iii) two open curves are equivalent if there exists 
an orientation preserving homeomorphism h of the plane such that 7' ([0,1]) = 
/i(7([0, 1])) and -f'{i) = h{-f{i)) for i = 0,l- 

Observe that, unlike closed curves, open curves are oriented from the initial 
point 7(0) to the final point 7(1). Moreover a unique closed curve is obtained from 
an open curve by connecting the final point to the initial one in counterclockwise 
direction around the curve. These definitions are illustrated by Figure H 
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Figure 2. A doodle and the associated rooted planar map. 

In order to study the families Ap and Ap and to obtain Conjecturenwe introduce 
a more general class of drawings, that we call doodles. For given positive integers p 
and k, let A^^p be the set of equivalence classes of (A:+l)-uples F = (70, 71, ... , 7fe) of 
curves drawn on the plane such that: (i) the curve 70 is an open curve of the plane: 
7o is a smooth mapping [0,1] M^, and 7o(0) and 7o(l) belong to the infinite 
component of \ (7o((0, 1)) U IJ^ ji{S^))- (ii) for i > I, each ji is a loop, that is 
a smooth map — > ; (iii) there are p points of intersection (including possibly 
self-intersections) of these curves, all of which are regular crossings; (iv) the union of 
the curves is connected, (v) two doodles F — (70, . . . , jk) and F' = {j'q, . . . , 7^) are 
equivalent if there exists an orientation preserving homeomorphism h of the plane 
such that 7^([0,1])uUtK51) = /i(7o([0, 1]) U U 7^(^')) and 7^(2:) = h{jo{^)), for 
a; = 0, 1. In other terms, a doodle is made of an open curve intersecting a set of 
loops, that are considered up to continuous deformations of the plane. An example 
of doodle is given in Figure [21 (left-hand side). 

2.2. Colored planar maps. An equivalent presentation of doodles is in terms of 
rooted planar maps |22[ \2'6\ . A planar map is a proper embedding of a connected 
graph in the plane considered up to homeomorphisms of the plane. It is ^-regular 
if all vertices have degree four. It is rooted if one root edge is marked on the infinite 
face and oriented in counterclockwise direction. Equivalently the root edge can 
be cut into an in- and an out-going half-edge (also called legs) in the infinite face. 
There is an immediate one-to-one correspondence between doodles with p crossings, 
and 4-regular planar maps with p vertices and two legs. This correspondence is 
illustrated on Figure |21 

We shall consider the number dAk^p of doodles with p crossings and 

k loops and more specifically we shall consider the asymptotic properties of ak,p as 
p (and possibly k) goes to infinity. It turns out to be convenient to introduce the 
generating function ap{n) as k varies: 

00 

(2) ap{n) y^Qfc,pn'' 

fc=0 

The requirement that a doodle is connected implies that it cannot contain more 
loops than crossings so that ap{n) is a polynomial in the (formal) variable^ n. 
For real valued n, ap{n) can be understood as a weighted summation over all 
doodles with p crossings, and, more specifically for n a positive integer, ap{n) can 
be interpreted as the number of colored doodles in which each loop has been drawn 
using a color taken from a set of n distinct colors. 

"'^The subsequent interpretation in terms of colored doodles and the strong tradition in the 
physics litterature are our (admittingly poor) excuses for the use of n to denote a formal variable. 
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On the one hand, m the special case fc = 0, ao,p — a-p gives by definition the 
number of open self- intersecting plane curves. Observe that a-p is also given by 
n — Q since ap(0) = ap^p. On the other hand, the generating functions of the 
ap{n) for other values of n, namely n = 1, 2, —2, have been computed exactly (see 
respectively EH 0^1 EH). We elaborate now on the case n = 1 since it will play a 
crucial role in what follows. 

The number ap(l) counts the number of doodles with p crossings irrespective of 
the number of loops k. In terms of maps, ap(l) is the number of rooted 4-regular 
planar maps with p vertices. The number of such planar maps is known [51 122j. 
from which one can compute the asymptotics: 

Observe in this case the power-law exponent —5/2, which is universal for rooted 
planar maps in the sense that it is observed for a variety of families of rooted planar 
maps (see (7j). As opposed to this, the exponential growth factor 12 is specific to 
the family of rooted 4-regular planar maps. 

There is a physical interpretation of the power-law behavior it is given 

by two-dimensional gravity. This explanation begs to be generalized to any n, and 
we shall explore such a possibility now. 



3. Two-dimensional quantum gravity and asymptotic combinatorics 

The purpose of this section is to give the rationale behind our conjectures. We 
place the discussion at a rather informal level that we hope achieves a double 
purpose: on the one hand it should give an intuition of the path leading to our 
conjectures to the reader with zero-knowledge in quantum field theory (QFT), and 
on the other hand it should convince the expert of this (quantum) field. For our 
defense, let us observe that filling in more details would require a complete course 
on QFT, with the result of not getting much closer to a mathematical proof. 

3.1. Prom planar maps to two-dimensional quantum gravity. The main 
idea of the physical interpretation of the numbers ap(l) is to consider planar maps 
as discretized random surfaces (with the topology of the sphere). As the number of 
vertices of the map grows large, the details of the discretization can be assimilated 
to the fluctuations of the metric on the sphere. To make this idea more precise, let 
us describe a way to associate a metric on the sphere to a given 4-regular map m: 
to each vertex of m associate a unit square and identify the sides of these squares 
according to the edges of m (arbitrary number of corners of squares get identified); 
the result is by construction a metric space with the topology of a sphere. Upon 
taking a 4-regular map uniformly at random in the set of map with p edges, a 
random metric sphere with area p is obtained. 

Now, physics tells us that the metric is the dynamical field of general relativity 
i.e. gravity, and that this type of fiuctuations in the metric are characteristic of a 
quantum theory. In our case it means that, as p becomes large, the discrete nature 
of the maps can be ignored and there exists a scaling limit, the properties of which 
are described by two-dimensional euclidian quantum gravity. In particular any pa- 
rameter of random planar maps that makes sense in the scaling should converge 
to its continuum analog. A fundamental parameter of this kind turns out to be 
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precisely the number of (unrooted) planar maps: it is expected to scale to the par- 
tition function Zg (A) of two-dimensional quantum gravity with spherical topology 
at fixed area A, through a relation of the form 

(4) -ap(l) ~ ZJA), with A proportional to p. 

p p^oc 

(Here the factor 1/p is due to the fact that the partition function does not takes 
the rooting into account.) The only thing we want to retain from Zg{A) is that 
the power law dependance of its large area asymptotic takes the form A~'^/'^, in 
accordance with Formula Q. (Trying to give here a precise description of the 
partition function Zg would carry us to far away, and anyway the arguments in this 
section are non rigourous.) 

In the case n ~ 1, this is the whole physical picture: a fluctuating but empty two- 
dimensional spacetime ~ there is no matter in it. What happens when n ^ 17 As 
already discussed, an appealing image is to consider that one must "decorate" the 
planar map by coloring each curve with n colors. Alternatively, the physicist's 
view would be to consider that we have put a statistical lattice model (of crossing 
loops) on a random lattice (the planar map). This perfectly fits with the previous 
interpretation of the planar map as a fluctuating two-dimensional spacetime: as 
we learn from physics, in the limit of large size, adding a statistical lattice model 
amounts to coupling matter to quantum gravity. Matter is described by a quantum 
field theory (QFT) living on the two-dimensional spacetime. The parameters of 
the lattice model that survive in the scaling limit are recovered in the critical (long 
distance) behavior of this QFT, which in turn is described by a conformal field 
theory (CFT). Then, provided we can find a CFT describing the lattice model 
corresponding to a given n 7^ 1, the analog of Relation Q holds with the partition 
function Zg_|_cFT(n)(^) of this CFT coupled to gravity: in the large size limit, 

(5) -ap{n) - Zg+cFT(«)(-4). 

In this picture, the only thing we need to know about the CFT that describes the 
scaling limit of our model is its central charge c, which roughly counts its number of 
degrees of freedom. Indeed, the study of CFT coupled to gravity was performed in 
IS El EI J resulting in the following fundamental prediction: the partition function 
Zg^CFT^A) of gravity dressed with matter has a power-law dependence on the area 
of the form A'^~^ where the critical exponent 7 depends only on the central charge 
of the underlying CFT via (for c < 1) 

V(l-c)(25-c) 

(6) 7= ^ . 

Returning to our asymptotic enumeration problem (not forgetting the extra fac- 
tor p which comes from the marked edge), we find: 

(7) ap{n) ^ e'^P+(7-2)logp + «; 

p — >oc 

where a, k are unspecified "non-universal" parameters, whereas the "universal" 
exponent 7 is given by Eq. © with the central charge c of the a priori unknown 
underlying CFT(n). The absence of matter, that is the case n = 1, corresponds to a 
CFT with central charge c = 0: one recovers 7 — 2 = —5/2 as expected from Eq. 
In general, all parameters in Eq. Q are functions of n; assuming furthermore that 
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their dependence on n is smooth in a neighborhood of n = 1, one can recover by 
Legendre transform of cr{n) the asymptotics of ak,p as k and p tend to infinity with 
the ratio k/p fixed. Observe finally that the knowledge of the CFT could give more 
informations. For instance, the irrelevant operators of the CFT control subleading 
corrections to Eq. |7J. 

3.2. The identification of two candidate models. We now come to the issue of 
the determination of the CFT for a arbitrary n. An observation made in |57], based 
on a matrix integral formulation, is that this CFT must have an 0{n) symmetry 
(for n positive integer - for generic n this symmetry becomes rather formal and 
cannot be realized as a unitary action of a compact group on the Hilbert space). 
There exists a well-known statistical model with 0{n) symmetry, a model of (dense) 
non-crossing loops 16 , whose continuum limit for |n| < 2 is described by a CFT 
with central charge 

(8) CI = 1-6(7^-1/7^)2 n = -2cos(^g), 0<g<l 

In |27| it was speculated that there is no phase transition between the model of 
non-crossing loops, which we call model I, and our model of crossing loops, and 
therefore the central charge is the same and given by Eq. ©. If this were the 
case, the study of irrelevant operators of this CFT would allow moreover to predict 
that subleading corrections to Eq. |(7J have power-law behavior for all |n| < 2 with 
exponents depending continuously on n. 

However, another scenario is possible. In ^^Ij it was suggested that 0{n) models, 
for n < 2, possess in general a low temperature phase with spontaneous symmetry 
breaking of the 0{n) symmetry into a subgroup 0{n — 1). This is a well-known 
mechanism^ in QFT (see e.g. |^ chapters 14, 30), which produces Goldstone bosons 
living on the homogeneous space 0{n)/0{n — 1) = 5"^^. In the low energy limit 
the bosons become free and the central charge is simply the dimension of the target 
space S*""^: 

(9) cii = n - 1 n <2 

For generic real n {n < 2) this is only meaningful in the sense of analytic contin- 
uation, but we assume it can be done and call it model II. This CFT possess a 
marginally irrelevant operator, leading to main corrections to leading behavior (UJ 
of logarithmic type i.e. in log logp ^^^^ 

It was furthermore argued in ^7] that the critical phase of model II is generic in 
the sense that the low-energy CFT is not destroyed by small perturbations - the 
most relevant 0(n)-invariant perturbation is the action itself, which corresponds 
to a marginally irrelevant operator for n < 2. On the contrary, the model I of 
non-crossing loops is unstable to perturbation by crossings; some numerical work 
on regular lattices (at n = 0) |H1 tends to suggest that it flows towards model II. 

Note that both Conjectures (O and supply the correct value c = for n = 1 
and c = 1 for the limiting case n — 2.^ Of course, in no way do we claim that 
these are the only possible scenarios which fit with known results - one might 



The Mermin- Wagner theorem, which forbids spontaneous symmetry breaking of a continuous 
symmetry in two dimensions, only applies to n integer greater or equal to 2. 

■^Actually, the two resulting c = 1 theories are not identical: the one from model I seems to be 
the wrong one, although this is a subtle point on which we do not dwell here. 
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have a plateau of non-critical behavior (c = 0) around n = 1, for instance; or two- 
dimensional quantum gravity universality arguments might not apply at all in some 
regions of n - but they seem the most likely candidates and therefore it is important 
to find a numerically accessible quantity which at least discriminates between the 
two conjectures. 



4. The general conjectures and a testable parameter 

The physical reinterpretation of doodles as a model on random planar lattices 
has led us to postulate that the weighted summation over doodles satisfies 

ap(n)~co(7i)T(n)P.p^(")-2^ 

with the critical exponent j{n) given in terms of the central charge c(n) by 



c(n)-l- V(l-c(n))(25-c(n)) 
7 n) = . 

Moreover we have presented two concurrent models which fix the value of c(n). 
Since negative values of n create additional technical difficulties (appearance of 
complex singularities in the generating function, cf |28p. we formulate the conjec- 
tures in the restricted range < n < 2: 

Conjecture 2 (Model 1). Colored doodles are in the universality class of dense 
non-crossing loops, so that for < n < 2, n ~ — 2cos(7rg), 1/2 < g < 1, 

c{n) = 1- 6(^-1/ 

Conjecture 3 (Model II). Colored doodles are in the universality class of models 
with spontaneously broken 0{n) symmetry, so that for < n < 2, 

c{n) = n — I. 

Observe that Conjecture implies Conjecture ^ for n = 0, while Conjecture |3 
would give c(0) = 1 - 6(\/2 - l/\/2)2 = -2 and 7(0) = -1. According to the 
discussion of the previous section, Conjecture |3| appears more convincing. In order 
to get a numerical confirmation, we look for a way to discriminate between the two. 

Since the model at n = 1 is much easier to manipulate, we look for such a 
quantity at n = 1. Of course the known value of the exponent 7(1) is a natural 
candidate but as already mentioned both conjectures agree on this: we propose 
instead the derivative of the exponent at n = 1, 

(10) 7' = T-, 

an \n=i 

The reason that it can easily be computed numerically is that it appears in the 
expansion of the average number of loops (/c)^ for a uniformly distributed random 
planar map with p vertices. Indeed one easily finds 

(11) ^^)p^T- logap(") = (T'p + i'\ogp + n' + o{l) 

dn\n=i P-+00 

Here we have assumed expansion ((TJ to be uniform with smoothly varying constants 
(T(ri), 7(n), K(n) in some neighborhood of n = 1, and written a' = 

k' = , Kin). 
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The conjectures El and 13 provide the following predictions for 7': 



(12) 



[3^^^0.413... in CFT I 
1^ = 0.3 in CFT II 



The quantity {k)^ is not known theoretically, so that we cannot immediately 
conclude in either direction. However it is possible to estimate it numerically using 
random sampling. 



In this section we present the algorithm we use to sample a random map from 
the uniform distribution on rooted 4-regular planar maps with p vertices. The 
problem of sampling random planar maps with various constraints under the uni- 
form distribution was first approached in mathematical physics using markov chain 
methods 12, 2 . However these methods require a large and unknown number 
of iterations, and only approximate the uniform distribution. Another approach 
was proposed based on the original recursive decompositions of Tutte 22 but has 
quadratic complexity PP, and is limited as well to p of order a few thousands. 

We use here a more efficient method that was proposed in [201 1^ along with 
a new derivation of Tutte's formulas. The algorithm, which we outline here in 
the case of 4-regular maps, requires only 0{p) operations to generate a map with 
p vertices and manipulates only integers bounded by 0{p). Moreover maps are 
sampled exactly from the uniform distribution. The only limitation thus lies in the 
space occupied by the generated map. In practice we were able to generate maps 
with up to 100 million vertices, with a generation speed of a million vertices per 
second. 

The algorithm relies on a correspondence between rooted 4-regular planar maps 
and a family of trees that we now define. A blossom tree is a planted plane tree 
such that 

• vertices of degree one are of two types; buds and leaves; 

• each inner vertex has degree four and is incident to exactly one bud; 

• the root is a leaf. 

An example of blossom tree is shown on Figure 13 By definition a blossom tree 
with p inner vertices has p + 2 leaves (including the root) and p buds. Observe that 
removing the buds of a blossom tree gives a planted complete binary tree with p 
inner vertices, and that conversely 3*^ blossom trees can be constructed out of given 
binary tree with p inner vertices. Since the number of binary trees with p inner 
vertices is well known to be the Catalan number (^^) , the number of blossom 



Let us define the closure of a blossom tree. An example is shown on Figure 13 
Buds and leaves of a blossom tree with p inner vertices form in the infinite face 
a cyclic sequence with p buds and p + 2 leaves. In this sequence each pair of 
consecutive bud and leaf (in counterclockwise order around the infinite face) are 
merged to form an edge enclosing a finite face containing no unmatched bud or leaf. 
Matched buds and leaves are eliminated from the sequence of buds and leaves in 
the infinite face and the matching process can be repeated until there is no more 
buds available. Two leaves then remain in the infinite face. 



5. Sampling random planar maps 



tree is seen to be 
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Figure 3. A blossom tree and its closure. Buds are represented 
by arrows. Dashed edges connect pairs of matched buds and leaves. 



Proposition 1 f[2(Jp. Closure defines a (p + 2)-to-2 correspondence between blos- 
som trees and rooted four-regular planar maps. In particular the number of rooted 
four-regular planar maps is 

2 3P f2p\ 
p + 2 p+1 \p J' 

This proposition implies that to generate a random map according to the uni- 
form distribution on rooted 4-regular planar maps with p vertices one can generate 
a blossom tree according to the uniform distribution on blossom tree and apply 
closure. A synopsis of the sampling algorithms is given below. An implementation 
is available on the web page of G.S. 

Random sampling of a rooted 4-regular maps with p vertices. 

Step 1. Generate a random complete binary tree Ti according to the uniform 
distribution on complete binary trees with p inner vertices. (This is done in linear 
time using e.g. prefix codes ^^.) 

Step 2. Convert Ti into a random blossom tree T2 from the uniform distribution 
on blossom trees with p inner vertices: independantly add a bud on each vertex in 
a uniformly chosen position among the three possibilities. 

Step 3. Use a stack (a.k.a. a last-in-first-out waiting line) to realise the closure 
of T2 in linear time: Perform a counterclockwise traversal of the infinite face until p 
buds and leaves have been matched; when a bud is met, put b into the stack; when 
a leaf £ is met and the stack is non empty, remove the last bud entered in the stack 
and match it with £. 

Step 4- Choose uniformly the root between the two remaining leaves. 



6. Simulation results 

The algorithm described in the previous section allows to generate random rooted 
4-regular planar maps with p vertices and two legs, with uniform probability. One 
can compute various quantities related to the map thus generated and then average 
over a sample of maps, as always in Montecarlo simulations. Here the main quantity 
of interest is the number of loops of the map. If we generate N maps of size p so that 
the ith map has fcp ^ loops, 1 < i < N, then ^ ^iLi kp,i has an expectation value 
of {k)p and a variance of jf{{k'^))p, where {k)p = _^logap(n) and {{k'^))p — 
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Table 1 . Numerical values ki of the average number of loops of 
maps with p = 2^ vertices. The error (standard deviation) on the 
last digit is given in parentheses. 
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25192.45(3) 


50381.35(5) 


100759.0(1) 


201514.3(2) 


403023.8(4) 


806043.2(7) 



-£^^^_^\ogap{ri) (the latter can of course itself be estimated as the expectation 

value of J2fLi fcp,j " NiN~i) According to Eq. 0, both {k)p and 

{{k'^))p are of order p for p large. However, we are interested in corrections to the 
leading behavior of {k)p which are of order \ogp, cf Eq. Hll|l . so that we need to 
keep the absolute error small. This implies that the size of the sample N should 
scale like p, or that the computation time grows quadratically as a function of p. 

In practice we have produced data for p = 2^ with £ < 24, the sample size 
being of the order of up to 10^. To ensure a good sampUng we used the "Mersenne 
twister" pseudo-random generator ^H], which is both fast and unbiased. The last 
few values of £ are only given to show where the statistical error begins to grow 
large due to limited memory and computation time. Let us call kg the numerical 
value found for (fc)p^2«- The results obtained are shown on tabled 

First, as a rough check of the asymptotic behavior, let us define ue = 2ki — ke+i- 
If expansion Hll() is correct, then ue must display an affine behavior as a function 
of £: Ue — {£— 1)7' log 2 + k' + 0(1/^). Indeed, as one can see on Fig.E] this is the 
case. 



-1 




Figure 4. The set of points ue = 2ki—ki+i as a function of logp = 
log 2 • £ with their error bars, as well as a proposed asymptote of 
slope 0.3. 

By comparison with the proposed asymptote it seems clear that 7' is close to 
0.3. To make this statement more precise, one can try to fit the set of the ke to 
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Table 2. Fits for the kg. is the mmimized weighted sum of 
squared errors. 
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0.04804366 


0.04804365 
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0.04804363 


0.04804365 
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n'p + 7' logp + k', where ^ ranges from I = £min to £ = (max = 24 and imm is 
varied. The resuhs are reported on Tab. |21 Unfortunately the confidence level 
remains fairly low until imin becomes so high that statistical error is huge, which 
tends to indicate strong corrections to the proposed fit. 

It is important to understand that if Conjecture 2 were true, then the corrections 
to asymptotic behavior would be power-law - starting with p~^^^ . This means that 
the procedure used in Tab.[21should converge quickly to the correct values of a' , 7', 
k' (to check this we have performed a similar analysis with a model of non-crossing 
loops on random planar maps and obtained fast convergence with high accuracy - 
2 digits on 7'). Here the range of values of 7' seems to be 0.29-0.34, far from the 
value predicted by Conjecture 2. It is therefore our view that the numerical data 
render Conjecture 2 extremely unlikely. 

On the other hand, the value 0.3 predicted by Conjecture 3 remains possible. 
The fluctuations observed even for very high p would be caused by the logarithmic 
corrections present in model II due to the marginally irrelevant operator, as men- 
tioned in Section 13.21 This operator is expected to induce a correction in 1/logp 
(which is in principle computable exactly using quantum field theory techniques, 
since it is universal; progress on this will be reported elsewhere), plus higher correc- 
tions, all of which remain significant in our range of data. This would also explain 
why it is so hard to extract useful information from the first few (exact) values of 
ap{n) given in [1(11 ITT] . 

In conclusion, and in view of the theoretical as well as numerical evidence, our 
belief is that Conjecture 3 is indeed correct. 

7. Variants and corollaries 

First observe that planar maps have in general no symmetries. More precisely 
the fraction of planar maps with p edges that have a non trivial automorphism 
group goes to zero exponentially fast under very mild assumption on the family 
considered . If this (very plausible) property holds then a typical closed curve 
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will be obtained by closing d different open curves, where d is the degree of the 
outer face. But the average degree of faces in any fixed 4-regular planar map is 
four. Thus the relation ap Aup. 

Second let us give a property illustrating the importance of the critical exponent 
7 as opposed to the actual value of r. A closed plane curve C is said to be a- 
separable, for < a < 1 a constant, if there exist two simple points x and y of 
C such that T \ {x,y} is not connected and both connected components contain 
at least p" crossings. The pair {x,y) is called a cut of C. In other terms, C is 
a-separable if it is obtained by gluing the cndpoints of two big enough open plane 
curves (up to homeomorphisms of the sphere). 

Corollary 1. Assume Conjecture^ is valid, and consider a uniform random closed 
plane curve Tp with p crossings. The probability that Tp is 1-separable decays at 
least like p'^ ^ p^^''''^ . More generally, z/ a > 1/(1 — 7) = (7 — •\/T3)/6 == 0.56, the 
probability that Tp is a-separable goes to zero as p goes to infinity. 

For comparison, 7 = —1/2 and 1/(1 — 7) = 2/3 for doodles, which are thus easier 
to separate. 

Indeed let us compute the expected number of incquivalent cuts of a closed plane 
curve with p crossings. When considered up to homeomorphisms of the sphere, close 
plane curves with a marked cut are in one-to-one correspondence with pairs of open 
plane curves. Hence, with a factor p for the choice of infinite face, 

13) p.x <cst.p.\ — = 0{pq< ). 

Qfp p^ ^ 

p'=q ^ p'=q 

In particular if (7 3> p^/^^^''') this expectation goes to zero as p goes to infinity. 

It is typical that in the computation of probabilistic quantities, like in Equa- 
tion p3|l . the exponential growth factors cancel, leading to behaviors that are driven 
by polynomial exponents. This explains the interest of in these critical exponents 
and gives probabilistic meaning to their apparent universality. As a final illustra- 
tion of this point let us present two variants of Conjecture ^ (Definitions of prime 
self-intersecting curves and alternating knots can be found in |10II14| .') 

Conjecture 4. The number a'p of closed prime self-intersecting curves withp cross- 
ings and the number a'p of prime alternating knots with p crossings lie in the same 
universality class as closed self-intersecting curves: there are constants t' , t" , c' , 
c" such that 

I I Ip ■7—2 II II Hp T— 3 

ap ^ c T ^ ■ p' , ap ^ c T ^ ■ p' , 
where 7 is given in Conjecture^ 

Observe that knot diagrams are naturally considered up to homeomorphisms of 
the sphere |10l I14| , while we have considered plane curves up to homeomorphisms 
of the plane. This explains the discrepancy of a factor p in Conjecture ^ for a'^, 
since one of the p-\-2 faces of a spherical diagram must be selected to puncture the 
sphere and put the diagram in the plane. 

8. Conclusion 

We have given arguments supporting Conjecturenfor the asymptotic number of 
plane curves with a large number of self-intersections, as well as the more general 
Conjecture|2| The numerical results provided in Section|Slsupport Conjecturenonly 
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indirectly since they are related to another specialization of Conjecture|21(derivative 
at 71 = 1 versus n — 0). However the alternative proposal is not compatible with 
either of these new numerical results (as is the case of Conjecture [21) or earlier ones. 

Our method to test the conjecture could be applied to other models like open 
curves with endpoints that are not constrained to stay in the infinite face, or the 
meanders studied by Di Francesco et al. 

Acknowledgements. P.Z.-J. would like to thank J. Jacobsen for pointing out 
Rcfs. ^ and jITl to him. 
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